function eq = find_agg_mit_ces(eq, param, glob, options)

    eq.C = (sum(eq.yf.^((param.sigma - 1)/param.sigma).*eq.L))^((param.sigma)/(param.sigma-1));

end